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INTRODUCTION 


Several  types  of  data  adaptive  detection  structures*  have  been  proposed  for  use  in 
environments  where  the  signal  and  noise  statistics  are  only  vaguely  known  and  where  insuf- 
ficient a priori  information  is  available  to  design  optimal  detectors.  Adaptive  mean  level 
detectors  have  been  implemented  to  provide  improved  performance  over  conventional 
nonadaptive  systems  when  the  noise  statistics  are  nonstationary  (ref.  1-7).  Other  forms  of 
adaptive  receivers  based  on  pattern  recognition  techniques,  decision-directed  feedback,  non- 
parametric  statistics,  and  sequential  processing  are  described  in  references  8-14.  Adaptive 
implementations  of  linear  prediction  filters  have  also  been  proposed  for  spectral  analysis, 
instantaneous  frequency  estimation,  data  compression,  adaptive  deconvolution,  system 
identification,  adaptive  noise  suppression,  and  speech  encoding  applications  (ref  1 5-29). 
The  computational  simplicity  of  adaptive  linear  prediction  (ALP)  implementations  based  on 
the  Widrow-Hoff  least-mean-squares  (LMS)  adaptive  algorithm  is  discussed  in  references 
16,18,20. 

Several  recent  papers  (ref.  17,21-25)  have  been  concerned  with  the  use  of  adaptive 
linear  prediction  filters  as  self-tuning  filters  for  the  detection  of  narrowband  signals  in  noise. 
In  the  present  study  an  L-weight  ALP  filter  is  used  to  estimate  the  frequencies  at  which  co- 
herent energy  is  present  in  the  input  data  sequence,  and  a K-point  spectrum  analyzer  is  used 
to  examine  either  the  transfer  function  of  the  adaptive  filter  weights  or  the  power  spectrum 
of  the  adaptive  filter  output,  as  illustrated  in  figure  1 . Appropriate  selection  of  the  pre- 
diction distance  (A  in  figure  1 ) on  the  basis  of  available  a priori  information  concerning  the 
expected  autocorrelation  of  the  noise  and  signal  components  allows  the  ALP  filter  to 
suppress  either  the  correlated  or  uncorrelated  noise  components  of  the  input  data.  This 
device,  pictured  in  figure  2,  has  been  called  an  “adaptive  line  enhancer”  (ALE)  due  to  the 
narrowband  signal  enhancement  which  results  from  the  broadband  noise  suppression  cap- 
ability of  a linear  prediction  filter  (ref.  1 7). 

The  purpose  of  this  paper  is  to  analyze  ALE  detection  performance  on  the  basis  of 
reciever  operating  characteristic  (ROC)  curves  for  inputs  consisting  of  a sinusoidal  signal  of 
known  frequency  but  unknown  phase  embedded  in  white  Gaussian  additive  noise  of  known 
power.  This  is  the  classical  signal-known-except-for-phase  (SKEP)  detection  problem  for 
which  the  optimal  detector  is  well  known  (14,30,31).  When  the  available  observation 
interval  is  long,  however  it  is  often  impractical  to  implement  the  optimal  detector;  thus, 
practical  spectrum  analyzers  typically  employ  incoherent  averaging  to  minimize  hardware 
requirements.  For  completeness,  the  performance  curves  obtained  for  the  adaptive  detector 
are  compared  with  those  of  both  the  optimal  detector  and  a commonly  used  nonadaptive 
detector  illustrated  in  figure  3. 

It  is  generally  possible  to  design  a nonadaptive  detector  which  is  superior  to  an 
adaptive  detector  in  a stationary  operating  environment.  However,  adaptive  detectors 
can  prove  advantageous  in  nonstationary  operating  environments  (ref.  1-7,24).  A goal  of 
this  paper  is  to  provide  baseline  performance  criteria  to  aid  in  selecting  appropriate  ALE 
time  constants,  filter  leng'  s,  and  other  design  parameters  for  narrowband  detection  appli- 
cations where  the  SKEP  election  problem  represents  an  appropriate  analytical  model  for 


•An  adaptive  detection  system  is  one  in  which  the  detector  structure  is  based  to  some  extent  on  the 
received  data  and  which  incorporates  a mechanism  to  vary  the  detector  structure  as  the  received  signal 
and  noise  characteristics  evolve  in  time. 
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Figure  3.  Processor  structures  for  coherent  DFT  and  incoherently 
averaged  detector,  (a)  Coherent  DFT.  (b)  Incoherently  averaged 
detector. 


the  input  signal  and  noise  conditions  in  which  the  detector  may  operate.  These  results  in  ' 

turn  allow  the  sensitivity  in  performance  to  be  determined  as  a function  of  the  adaptive  j 

filter  parameters.  j 

The  Widrow-Hoff  LMS  algorithm  provides  a computationally  efficient  method  for  I 

updating  the  ALP  filter  coefficients  to  track  temporal  variations  in  the  signal  and  noise  * 

statistics  (ref.  1 6-20).  All  performance  analyses  of  the  ALE  in  this  paper  are  based  on  an 
ALP  filter  which  is  updated  by  this  algorithm  (figure  2). 

In  suceeding  sections,  the  sensitivity  of  the  ALE  detector  to  the  feedback  para- 
meter p is  investigated  and  found  to  be  relatively  small.  It  is  further  shown  that,  for  L = K, 
the  detection  performance  of  the  two  adaptive  detectors  in  figure  1 compares  favorably 
with  that  of  N/K  independent  incoherent  averages  of  a K-point  DFT  (figure  3),  provided  N, 
the  length  of  the  data  sequence  is  of  the  order  of  the  ALE  convergence  time.  Finally,  a 
study  is  made  of  the  effects  of  varying  K. 
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ADAPTIVE  WEIGHT  VECTOR  MODEL 


As  illustrated  in  figure  1 , the  ALE  consists  of  an  L-weight  linear  prediction  filter  in 
which  the  filter  coefficients  are  updated  at  the  input  sampling  rate.  The  adaptive  filter  out- 
put r(k)  is  defined  by 


L-1 

r(k)=  2 Wj(k)x(k-j-A)  , (1) 

j=0 

where  A is  the  prediction  distance  of  the  filter.  The  output  r(k)  is  then  subtracted  from 
x(k)  to  form  an  error  sequence  e(k).  The  error  sequence  is  fed  back  to  adjust  the  filter 
weights  according  to  the  Widrow-Hoff  LMS  algorithm  (ref.  30), 

Wjfk  + 1)  = Wj(k)  + 2p  e(k)  x(k -i  - A),  i = 0,1, . . .,  L -1  , (2) 

where  p is  a constant.  Tlie  ALE  weight  vector  w(k  + 1)  at  time  k + 1 may  be  written  from 
(2)  as 


w(k  + 1 ) = w(k)  + 2/i  x(k)  [x(k)  - x^(k)  w(k)]  (3) 

where 

x(k)  = [x(k  - A),  x(k  - A - 1),  . . .,  x(k  - A - L + 1)]^ 

and 

w(k)  = (wQ(k),  W|(k), . . .,  WL_j(k)]^  . 

Taking  the  expectation  of  both  sides  of  (3),  we  have 

E{w(k  + 1)}  = E{w(k)}  + 2p[E{x(k)  x(k)} -E{x(k)  x^(k)  w(k)}]  . (4) 

Since  E[w(k  + 1)]  depends  on  all  the  past  data  x(k),  x(k  - 1),  . . .,  computation  of  the  con- 
vergence of  w(k)  in  the  mean  becomes  intractable  unless  some  simplifying  assumptions  are 
made.  For  small  p and  low  inpul  signal-to-noise  ratio  (SNR),  there  is  very  little  statistical 
dependence  between  the  x(k)  and  w(k)j  i.e., 

E[x'r(k)  w(k)l  s E[(xT'(k)]  E(w(k)J  , 

and  (5) 

E {x(k)  x^(k)  w(k)}  = E{x (k)  xT(k)}  E {w(k)}  . 

Widrow  (ref.  30,1 7)  has  shown  that,  for  uncorrelated  data,  E{w(k)}  converges  from  an 
arbitrary  initial  value  to  w*  as  k -*•  «>.  provided  0 < p < where  is  the  largest 

eigenvalue  of  the  data  autocorrelation  matrix  and  w*  is  the  optimal  Wiener  filter  weight 
vector.  The  term  w*  can  be  obtained  from  the  solution  to  the  discrete  Wiener-Hopf  matrix 
equation. 


where  R is  the  L X L autocorrelation  matrix  of  the  input  sequence  x(k)  and  P is  the  L-element 
cross-correlation  vector  with  elements  Pj  = E {x(k)  x(k  + j + A - 1)}.  The  convergence  of  the 
mean  weight  vector  in  the  more  general  cases  of  correlated  data  and  for  nonstationary  inputs 
has  been  considered  for  both  ideal  filters  and  for  cases  in  whch  implementation  errors  are 
present.  It  is  shown  that  the  convergence  properties  of  the  filters  are  robust  under  a wide 
range  of  conditions,  providing  that  the  inputs  to  the  filter  satisfy  specified  nondegeneracy 
conditions  (ref.  33-34).  Experimental  measurements  of  the  mean  weight  vector  of  a hardware 
implementation  of  the  ALE  also  confirm  that  E{w(k)}  -*■  w*  for  increasing  k for  sinusoidal 
inputs  in  uncorreclated  noise  for  many  cases  of  interest  (ref.  23). 

Since  the  LMS  algorithm  only  approximates  the  true  LMS  error  gradient,  the  weight 
vector  coefficients  are  obtained  from  noisy  estimates  of  the  optimal  Wiener  filter  solution. 
This  gradient  estimation  noise,  often  called  “misadjustment  noise”  (ref.  17),  is  given  by 

Vj(k)  = Wj(k)  - w*  . (7) 

Widrow  et  al.  (ref.  17)  have  shown  that  the  variance  of  Vj  may  be  made  arbitrarily  small  by 
decreasing  the  feedback  parameter  //.  Decreasing  ti  also  increases  the  convergence  time  of 
the  LMS  algorithm,  however,  and  a compromise  between  the  response  time  and  the  misad- 
justment noise  power  must  be  made  in  practical  applications. 

It  is  generally  assumed  (ref.  1 7)  that  the  misadjustment  noise  has  a Gaussian  distribu- 
tion. It  then  follows  from  (5)  that 

E[Vj(k)]  =0  , 

(8) 

E[Vj(k)Vi(k)]  =M?„in6(i-j)  , 

where  6(')  is  the  Kronecker  delta  operator  and  is  the  minimum  mean  square  error 
power.  Analytical  solutions  for  W:  for  multiple  sinusoids  in  uncorrelated  noise  have  been 
derived  (ref.  23).  When 

x(k)  = A sin  (coQk  + &)  + n(k)  , (9) 

where  n(k)  is  a zero-mean  white  Gaussian  noise  sequence  of  variance  v~,  and  wq  is  known 
a priori,  the  sampling  frequency  can  be  adjusted  so  that 

COqL 

— — = integer  , (10) 

2n 

where  cog  is  normalized  relative  to  the  sampling  frequency  of  lir  radians.  In  this  case  the 
R-matrix  is  circulant,  and  it  can  readily  be  shown  (ref.  21,23)  that 

w*  =-^  cos  [cjQ(j  + A)l,  0 <j  < L- 1 , (11) 


where 


^ (L/2)  SNR 
1 + (L/2)  SNR  ’ 
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and  the  input  SNR  is  defined  by 
SNR  = a2/2u2  . 


(13) 


The  amplitude  of  w.*  depends  nonlinearly  on  the  input  SNR,  since  the  transversal  filter  can- 
not discriminate  between  the  sinusoid  and  the  broadband  noise  within  the  filter  passband. 

It  has  been  shown  (ref  23)  that  (11)  also  provides  an  excellent  approximation  to  w*  for 
(jt/L)  < cuo<  [L  - 1/(L)1  for  noncirculant  autocorrelation  matrices. 

The  speed  of  adaptation  of  the  LMS  algorithm  is  controlled  by  the  dominant  eigen- 
values of  the  R-matrix  (ref  21,34).  For  the  case  of  a single  sinusoid  in  uncorrelated  noise, 
the  optimal  weight  vector  is  orthogonal  to  all  eigenvectors  of  R except  the  conjugate  pair 
corresponding  to  frequency  u>q.  The  eigenvalues  associated  with  these  eigenvectors  dominate 
the  convergence  time  and  are  both  approximately  given  (ref.  21)  by 


X = v^  + (A^L/4)  . 


(14) 


When  the  filter  is  filled  with  data  and  all  weights  are  zeroed  at  k = 0,  the  mean  value  Wj(k) 


of  the  jth  weight  is  updated  according  to  the  relation 
Wj(k)=  [1 -(1 -2mX)'^]  w?,  k = 0,l,... 
For  2juX  « 1 and  k » 1,  (15)  is  approximated  by 
Wj(k)  = [ 1 - w* 

= ^^^cos[wo(j  + A)],0<j<L  - 1 


(15) 


where 


a(k)=  (1  -e  '"la* 


(16) 


(17) 


In  (17),  Tjjj  = 1/2/iX,  the  mean  adaptation  time  constant. 

It  was  shown  (ref  35)  that  the  power  spectral  density  of  the  ALE  output  for  an 
input  consisting  of  a sinusoid  in  white  noise  contains  three  distinct  spectral  components:  a 
sinusoidal  term  due  to  the  filtering  of  the  input  sinusoid,  a white  noise  term  due  to  the 
misadjustment  noise,  and  a narrowband  noise  term  due  to  the  white  noise  passed  through 
the  converged  Wiener  filter.  The  relative  magnitudes  of  these  components  were  derived  in 
reference  35.  It  was  further  shown  that 


2n^.,2 


•min 


(18) 


We  now  derive  the  detection  performance  of  the  two  distinct  implementations  of 
the  ALE  illustrated  in  figure  1 . ROC  curves  are  obtained  under  the  condition  that  the 
adaptive  filter  weights  are  initialized  to  zero  at  k = 0,  and  that  N data  points  are  processed 
prior  to  forming  a detection  statistic.  The  L-weight  adaptive  filter  and  the  K-point  DFT  are 
restricted  so  that  K and  L are  much  less  than  N.  This  restriction  eliminates  the  initial  trans- 
ient behavior  as  the  filter  fills  with  data  and  permits  the  mean  weight  approximation  (16). 

In  all  cases  A = 1 . 
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DETECTION  STATISTICS  OF  ALE  OUTPUT  DETECTOR 


The  ALE  output  detector  is  shown  in  figure  1(b).  Detection  is  based  on  a K-point 
DFT  of  the  final  K output  values  from  the  L-length  filter  using  the  classical  Neyman-Pearson 
hypothesis  test.  Define  the  two  hypotheses,  Hq  and  Hj,  by  Hq  : x(k)  = n(k)  and  Hj : x(k) 

= s{k  I d)  + n(k),  where  s(k  | 0)  is  the  sinusoidal  signal  of  known  amplitude  and  frequency 
but  unknown  phase  and  n(k)  is  a sample  of  the  Gaussian  uncorrelated  noise  sequence  with 
zero  mean  and  variance  i.e., 


s(k  1 0)  = A sin  (wQk  + 0)  , 

(19) 

E[n(k)n(p)]  =u^6(k-p)  , 

(20) 

where  A,  coQ.and  are  known  a priori  and  6 is  unknown. 

Since  s(k  1 d)  is  zero-mean,  the  variable  x(k)  is  a zero-mean  Gaussian  random  variable 
for  all  values  of  k.  Thus  the  filter  output  r(k  | d)  conditioned  on  the  unknown  initial  phase 
defined  by  (1)  is  a linear  combination  of  products  of  independent  Gaussian  variables.  It  will 
be  assumed  that  for  large  L,  r(k  1 0)  is  approximately  Gaussian.  The  output  spectrum 
obtained  from  a K-point  DFT  of  r(k|  6)  evaluated  at  cjq  is  given  by  RIcoq)  = u + jv,  where 
j =\/^  and 

K-1 

u = ^ r(k  I 0 ) cos  coQk  , 
k=0 

K-1 

v=  ^ r(k  I 0)  sin  wok  . (21) 

k=0 


The  detection  statistic  q is  then  defined  as 

q = u^  + v^  . (22) 

The  joint  probability  density  function  (PDF)  for  (u,v)  must  therefore  be  specified  in  order 
to  define  the  PDF  of  the  detection  statistic  q.  It  is  shown  in  the  appendix  that  u and  v are 
essentially  statistically  independent  for  the  given  signal  and  noise  models  above.  Since 
r(k  I 0)  is  assumed  Gaussian,  u and  v are  also  Gaussian  and  the  PDF  p(u,v)  may  be  written  as 

p(u.v)M2.v,J-'  I . (23, 

where  u and  v are  the  component  means  and  and  their  variances. 

Expressions  for  u,  v,  a^,  and  are  derived  in  the  appendix,  and  appear  in  equations 
(A.9),  (A.  10),  (A.52),  and  (A. 53),  with  the  attendent  restriction  that  (tt/M)  « wq  « tt 
- (jt/M),  where  M = min  (K,L).  (See  references  36  and  37  for  additional  detail.)  Provided 
that  N > Tj^  and  that  (L,K)  « N,  a*  in  (A.9),  (A.  10),  and  (A.52)  can  be  approximated  by 
a(k),  as  defined  by  (17),  to  describe  the  output  statistics  just  prior  to  convergence.  Thus, 
under  H j : 

u = A sin  0;  V = A cos  0 , (24) 


I' 

(i 
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r 


where 


A = a(k)  AK/2. 

2 

In  addition,  the  variance  under  Hj,  Oj,  can  be  expressed  as 

2 2.  _2  _ ^^min*^^^^  . ^^min-^  j.  2a^(k)i>^D 

<^l-«u-^v-  2 8 t2 


where 


K^L 

4 

12 

L^K 

l3 

4 

12 

The  first  two  components  of  are  independent  of  the  time  index  because  of  the  assumption 
that  the  weight  misadjustment  noise  is  independent  of  k in  (8).  The  first  term  in  (26)  is  due 
to  the  broadband  weight  misadjustment  noise  produced  by  the  input  white  noise.  The  second 
component  of  arises  from  a narrowband  noise  process  centered  at  frequency  ccq  due  to 
the  weight  misadjustment  caused  by  the  signal.  This  narrowband  noise  process  may  be 
intrepreted  as  a random-amplitude  sinusoid  with  frequency  coq.  The  third  component  of 
the  variance  arises  from  the  white  noise  passed  by  the  adaptive  filter  and  increases  in  mag- 
nitude with  increasing  k.  Under  Hq,  a(k)  = 0,  and  the  component  means  and  variances 
defined  by  (24-26)  reduce  to  iT  = v = 0 and 


^^min  ^ 


The  derivation  of  the  probability  density  function  for  the  detection  variable  q defined 
by  (22)  is  straightforward  (ref.  14,30).  It  is  readily  shown  that,  under  Hq,  pQ(q)  is  the  gamma 
density  function  with  two  degrees  of  freedom  given  by 


P0(q)=  exp  r-q/2a^]  , q>0  . 
2a5 

Under  Hj,  q has  a density  function  pj(q)  given  by 


q>0  , 


where  Iq(x)  is  the  modified  Bessel  function  of  the  first  kind. 

According  to  the  Neyman-Pearson  criteria,  the  probability  of  false  alarm  Pp^  deter- 
mines the  threshold  under  Hq: 


10 


PfA=  / • 


(30) 


Since  Pp^  and  po(q)  are  known,  (3)  completely  specifies  the  threshold  Pq.  When  the  signal 
is  present,  the  probability  of  detection  Pd  is  given  by 


Pd^  f Pi(q)dq  ■ 

<i=Po 


(31) 


Through  a change  of  variables  it  may  be  shown  (ref.  30)  that  equivalent  forms  of 
these  integrals  are 


PfA=  f X5(p)dp 

P=To 


(32) 


and 


D 


= f 4 


(p  I a)  dp  , 


(33) 


P=T, 


where  X2(P)  is  the  chi-square  PDF  with  2 degrees  of  freedom  and  X2'  (P  I a)  is  the  non- 
centrality parameter  a,  with 


oi  = A^/a^  = 


K^a^fk)  SNR 


Lfii 


min 


r ^ K-SNR  ^ 4a2(k)D] 

L 2 . l3  J 


(34) 


Specification  of  Pp^  in  (30)  determines  the  threshold  Tq.  This  threshold  is  then  scaled  by 
the  ratio  of  the  variances  under  Hq  and  Hj  to  obtain  the  threshold  Tj 


Ti  - (aQ/aj)2  Tq  - 


+ -i^SNR-H  4a^(k)D  1 ~*T 

2 uf  . L%  ° 

^smin^  ^ J 


(35) 


The  probability  integrals  (30,  31)  may  be  evaluated  by  a nonnumerical  method  developed  by 
Urkowitz  (ref.  38)  or  by  standard  techniques  for  evaluating  Marcum  Q-functions  or  incom- 
plete Toronto  functions. 


r 


DETECTION  STATISTICS  OF  THE 
ALE  WEIGHT  VECTOR  IMPLEMENTATION 


Reeves  (ref.  39)  has  derived  ROC  curves  for  the  SKEP  detection  problem  using  the 
ALE  weight  vector  implementation  illustrated  in  figure  1(a).  Reeves  bases  detection  on  a 
single  L-point  ALE  weight  vector  after  the  ALE  has  processed  N data  points,  where  L « N. 
The  detection  statistic  is  defined  as 


9 9 

q = u,,,  + v,.„ 
W w 


(36) 


where 


%- 


L-1 

^ Wg(k)  cos 
C=0 


and 


'^w- 


L-1 

^ wg(k)  sin  coqC 

e=o 


(37) 


It  is  assumed  (ref.  26)  that  w£(k)  is  given  by  (16)  and  that  the  weight  variance  (equal  under 
Ho  and  Hj)  is  given  by  (8).  Therefore  PpA  and  Pd  can  be  obtained  from  integrals  of  the 
form  (ref.  26): 


oo 

PfA=  f xl(z)dz 


(38) 


Z=Tr 


and 


>0=  / 

z=T, 


xj  (z 


oi^)  dz 


(39) 


0 


where  = 2?(k)/a,^,,  and  is  the  weight  variance  defined  by  (8);  i.e 

W a W W 


Of..  = Mfr 


2 

..  „ „ „ w '*'min‘ 

(Note  that  the  lower  limits  of  the  Pq  and  Pp^  integrals  are  equal  in  this  case.)  These 

integrals  (38,39)  can  also  be  evaluated  using  the  nomograms  in  reference  38. 


12 


OPTIMAL  DETECTION 


L. 


jl 

The  optimal  likelihood  ratio  detector  for  a sinusoid  of  unknown  amplitude  and  phase 
is  shown  in  figure  3 (ref.  30).  Since  the  frequency  of  the  sinusoid  coq  is  a known  priori,  this  ; 

optimal  detection  structure  can  be  implemented  using  a DFT  which  transforms  the  input 
sampled  data  sequence  to  ensure  that  tJQ  is  one  of  the  DFT  frequencies.  The  real  and 
imaginary  parts  of  the  DFT  are  squared  and  summed  to  form  the  test  statistic.  The  test 
statistic  is  formed  from  a coherent  summation  of  the  squared  output  over  the  entire  observ- 
ation interval.  For  the  Gaussian  white  input  noise  sequence  defined  by  ( 1 ),  the  resulting 
Pp^  and  Pq  are 

OO 

PpA  = J X2(z)  dz  (40) 

z=/3o 

and 

OO 

Pd=  / X2’(z|NSNR)dz  , (41) 

where  N is  the  total  number  of  points  in  the  observation  interval. 


DETECTION  STATISTICS  FOR 
INCOHERENTLY  AVERAGED  DFT  PROCESSOR 

In  applications  where  the  total  time  duration  of  the  signal  N is  unknown  a priori,  or 
in  cases  where  the  input  time-bandwidth  product  is  large,  it  may  not  be  practical  to  match 
the  length  of  the  coherent  summation  to  N.  In  such  cases,  the  suboptimal  structure  shown 
in  figure  3(b)  is  often  used.  In  this  case  the  N-length  input  sequence  is  divided  into  M blocks 
of  K points  each,  and  a coherent  K-point  DFT  is  formed  for  each  block.  The  block  DFTs 
are  then  averaged,  forming  the  detection  statistic  u (normalized  to  unity  variance),  which  is 
the  magnitude  square  of  the  cjq  component.  Under  Hq  the  detection  statistic  u is  a chi- 
squared  density  function  with  2M  degrees  of  freedom  (ref.  30,38).  Thus 

OO 

PpA=  f • ^"^2) 

Under  Hj  the  detection  statistic  u is  a noncentral  chi-squared  density  function  with  2M 
degrees  of  freedom  and  noncentrality  parameter  Oj  = N(SNR)  (ref.  38).  Since  the  statistical 
variance  on  u is  not  dependent  upon  signal  presence  or  absence,  the  threshold  under  Hj  is 
also  equal  to  Pq.  Thus 

OO 

Pp)  = J'  X2M  I “l^u  . (43) 

M=|3o 
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The  quantities  |3q  and  Pq  may  be  evaluated  using  the  nomograms  presented  in  reference  38. 
Note  that  the  performance  of  the  quadratic  detector  for  this  case,  as  shown  in  figure  6 below, 
is  closely  approximated  by  the  ROC  curves  developed  by  Robertson  (ref.  40).  An  extensive 
discussion  of  the  performance  of  this  class  of  detectors  is  presented  in  reference  30. 


PERFORMANCE  ANALYSIS 


The  detection  performance  of  the  adaptive  implementations  and  the  conventional 
DFT  processor  of  the  previous  section  will  be  compared  on  the  basis  of  ROC  analysis.  A 
value  of  Ppy^  is  specified  and  held  constant  for  a set  f'f  hypotheses  tests,  while  the  input 
signal  strength  is  varied.  The  results  are  sets  of  curves  which  give  Pq  as  a function  of  input 
SNR  for  different  processor  configurations. 

For  the  ALE  detectors  there  is  a fundamental  trade-off  in  the  selection  of  p be- 
tween the  midadjustment  noise  (8)  and  the  convergence  time  (17).  In  the  output  detector, 
for  a given  Pp^,  the  probability  of  detection  is  an  increasing  function  of  a and  a decreas- 
ing function  of  oq/uj  (26,27).  Furthermore,  from  (34,35),  it  can  be  seen  that  a and  oq/O] 
are  monotonic  functions  (increasing  and  decreasing  respectively)  of  a“(N)/p,  where  N is 
the  length  of  the  data  sequence.  Similarly,  Reeves  (ref.  39)  has  shown  that  the  detection 
performance  of  the  weight  implementation  is  a function  of  a^(N)/p.  Thus,  maximizing 
detectability  at  time  N with  respect  to  p is  equivalent  to  maximizing  a^(N)/p  and  yields  the 
following  relationsjiip  between  p^p  and  N for  both  adaptive  implementations: 


1.25643 


op 


N 2u“ 


1 + 


SNR 


N»L. 


(44) 


It  follows  from  ( 1 7)  that  optimal  detection  is  obtained  at  a(k)  = 0.7 1 5a*,  indicating  that 
the  ALE  weight  vector  has  reached  71 .5  percent  of  its  steady-state  value. 

Detection  performance  is  not  critically  dependent  on  the  choice  of  p when  p is  close 
to  Pgp.  This  is  illustrated  in  figure  4,  where  the  SNR  at  which  Pd  = 0.5  is  plotted  as  a 
function  of  p for  Pp^  = 10~^. 

Figure  5 displays  Pj^  versus  input  SNR  for  the  case  K = L = 1024  and  an  input  data 
sequence  of  N = lOOL  points.  The  ROC  curves  for  the  ALE  implementation  are  compared 
with  those  of  a coherent  DFT  of  length  N and  an  incoherently  averaged  DFT  (lAD)  proces- 
sor which  averages  the  power  spectra  of  1 00  K-point  DFTs.  The  performance  of  the  two 
adaptive  detectors  is  nearly  equivalent  (within  0.5  dB)  to  that  of  the  lAD  process  in  the 
vicinity  of  Pp  = 0.50. 

Figure  6 shows  the  effect  of  the  DFT  length  on  detection.  The  primary  difference 
from  figure  5 is  a rise  in  the  detectable  SNR  levels  reflecting  the  loss  in  coherent  gain, 

10  log  K|/K2,  caused  by  the  decreased  resolution.  On  the  other  hand,  changing  the  ratio 
N/L  affects  the  incoherent  gain,  producing  an  approximate  loss  of  5 log  ((N  j L2)/(N2L])1 
(compare  figures  5 and  7). 

Note  that  the  value  of  p was  optimized  in  the  above  cases.  The  relative  behavior  of 
the  lAD  and  the  two  adaptive  detector  implementations  also  remained  similar  when  the  DFT 
and  ALE  lengths  were  shortened  to  K = L = 128. 
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Figure  4.  Sensitivity  in  detection  performance  to  choice  of 
adaptive  time  constant  for  ALE  weight  and  output  detectors; 
SNR  versus  ju  for  Pq  = 50  percent,  =■  10"^,  K = L and 
N = lOOL.  Under  these  parameters  the  ALE  weights  and 
ALE  output  differed  by  less  than  0.2  dB,  and  thus  a single 
curve  appears  for  both  detectors. 


Figure  5.  ROC  curves  for  PpA  = lO"'^,  N = 102400,  K = L = 
1024,  f/i  min  = 4. 1 X J0"h.  ^ jj  optimal  at  SNR  = 30.1  dB 
(Pjj  ^ 50  perccni). 
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Figure  6.  KOC  curves  for  Pp^  = 10"^,  N = 12800,  K = L = 
128,  /i^  rnin=  3.3  X 10"^;  is  optimal  at  -21.2  dB  (Pjj  SO 
percent).  The  ALE  weights  and  lAD  coincided  within  the 
accuracy  of  the  figure. 


Figure  7.  ROC'  curves  for  1^^  = 10"'^,  N = 20480,  K = L = 
1024,  = 1.54  X 10-5;  is  optimal  at  SNR  = 27.1  dB 

(Pjj  3 .SO  percent). 


Figure  8 illustrates  the  effects  of  changing  K/L  for  the  ALE  output  device.  ALE 
detection  performance  is  a nonlinear  function  of  K/L  as  indicated  by  (34)  and  (35).  In 
figure  8,  L is  fixed  at  1024,  and  K is  variable  from  64  to  16,192.  The  parameter  p is  op- 
timized at  each  SNR  using  (44).  A variation  of  between  3 log  (K/L)  and  5 log  (K/L)  is 
indicated  as  K/L  is  varied  by  factors  of  4. 


Figure  8.  ROC  curves  for  ALE  output  detector  for  fixed  fil- 
ter length  L = 1 024  and  varying  DFT  length  K with  = 
10*4,  = 4. 1 X 10*6,  and  N = 102400;  p is  optimal  at 

SNR  = 30. 1 dB. 


SUMMARY 

This  paper  has  examined  the  detection  properties  of  two  adaptive  detectors  for  the 
SKEP  probelm.  ROC  curves  were  computed  for  the  two  adaptive  implementations  (an 
L-point  ALE  followed  by  a K-point  DFT  of  the  filter  weight  vector  or  output)  and  compared 
to  those  of  an  optimal  DFT  of  length  N and  a K-point  DFT  with  coherent  averaging.  The 
effects  of  the  parameters  K and  L and  of  the  data  sequence  length  N were  discussed.  It 
was  found  that,  when  K = L,  the  two  adaptive  detector  implementations  performed  similarly 
to  an  incoherently  averaged  DFT,  provided  that  the  feedback  constant  p was  optimized  for 
a given  N,  L,  and  SNR.  The  sensitivity  of  detection  performance  to  the  choice  of  the  adap- 
tive time  constant  for  a given  observation  interval  was  discussed.  Finally  the  effect  of  a 
resolution  mismatch  between  the  ALE  and  the  DFT  was  investigated. 

The  results  obtained  in  the  paper  are  based  on  the  assumptions  that  the  ALE  output 
noise  statistics  are  Gaussian  and  that  the  magnitude-square  transform  of  the  ALE  output  has 
a chi-square  density  with  two  degrees  of  freedom.  Future  work  will  concentrate  on  defining 
these  densities  experimentally  and  obtaining  experimental  verification  of  the  analytical  ROC 
curves  presented  here.  It  will  also  include  the  effects  of  post  detection  integration. 
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APPENDIX 


The  purpose  of  this  appendix  is  to  calculate  the  means  and  variances  of  the  two 
variables  u and  v appearing  in  equation  (23). 

Two  expressions, 

M-1 

Gi^  = Gj  (CO.  b)=  cos(2wk  + b)  (A.l) 

k=0 


and 

M-1 

G!^.  = -L—  fcj,  b + -^^=  kcos(2cok  + b)  , 

” 2 V 2 ) ^ 

will  occur  frequently  throughout  the  calculations.  These  expressions  may  be  evaluated 
(using  standard  identies)  through 

G smJlco^in(^(M_i)+b)  , (A.2) 

“ sin  CO 


and  are  bounded  (ref.  33)  by 
I G^^  I « M 

I G'^^  1 « m2  (A.3) 

under  the  condition 


— « CO  « IT  - — . (A.4) 

M M 

Throughout  the  remainder  of  the  appendix  (unless  otherwise  stated)  we  shall  assume  that 
(A.4)  is  satisfied  where 

M = min(L.K)  . (A.5) 

In  addition,  we  will  assume  that  the  delay  is  1 : 

A = 1 (A.6) 

Two  other  identities  will  be  useful: 

cos  co(£  - k)  sin  coC  cos  cok  = — [sin  2co(C  - k)  + sin  2co£  + sin  2cok]  (A.7a) 


and 

cos  co(C  - k)  cos  coC  cos  cok  = — ( 1 + cos  2co(6  - k)  + cos  2coC  + cos  2cok]  , (A.7b) 

4 

which  may  easily  be  verified  by  substituting  exponentials;  e.g.,  cos  coC  = y (e'^*^®  + e~‘^®),  etc. 


JL  la 
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CALCULATION  OF  MEANS 


K-1  L-1 

E(u)=  E[Wg(k)X(k-6- Dcoscokl 

k=0  JP=0 

= A 2 Wj  sin  (a)(k  - C - 1)  + 0)  cos  wk 
k,C 


2a*A 


sin  (cj(k  - C - 1)  + 0)  cos  o;(K  +1)  cos  cok 
k,8 


= - 2 sin  [(co(k  - C - I ) cos  0 + cos  co(k  - C - 1 ) sin  0 ] 

k,e 

[cos  (cc(6  + 1 ))  cos  cok]  . 

Use  of  (A.7)  and  the  identity  sin  (con  + 0)  = sin  con  cos  0 + cos  con  sin  0 yields 

E(u)  = [sin  0 + sin  (2co(k  - 8 - 1)  + 0)  + sin  (2cok  + 0)  + sin  (2co8  + 0] . 

“ M (A.8) 


Then,  from  (A.3), 


L-1 


E(u)  = ^ (K  sin  0 + terms  « K)  s 

8=0 


a*AKL 

2L 


sin  0 . 


Similarly 


E(v)s 


a*AK 


cos  0 


(A.9) 


(A.  10) 


More  exact  expressions  may  be  found  in  reference  37. 


CALCULATION  OF  aj  AND  oj 

The  variance  of  u is  given  by  (23); 
ol  = E(u2)  - E(u)2 
K-1 

= 2 7^(8  - k)  cos  cok  cos  co8  , (A.ll) 

8,k=0 
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where 


f 


I 


7j.(m)  = E{(r(k  I 0)  r(lc  + m I 0)}  - E(r(k  1 0)1  Elr(k  + m I 0)1  . (A.12) 

It  follows,  from  the  approximate  independence  of  Wj  and  Xj  (5),  that  for  A = 1 

E(r(k  I 0))  = E I ^ wg(k)  x(k  - C) 

' S=0 

L-1 

= 2 E[(wg(k)l  E[Asin((k- 

e 

L-l 

= ^ wg  A sin  (aj(k  - C - 1)  + 

C 

Now  let 

0(m)  = E[r(k|0)r(k  + m|0)l 
L-l  L-l 

= E ^ [wj(k)Wj(k  + m)][Asincj(k-j- l)  + 0]  + [n(k-j- A)1 
j=0  i=0 

• [A  sin  to(k  + m - i - 1)  + 01  + [n(k  + m - i - 1)1  . (A.14) 

Again,  since  the  weight  vectors  wj  and  the  white  noise  input  n are  essentially  uncorrelated, 
and  since  E{n)  = 0,  most  of  the  terms  in  (A.  14)  vanish.  Two  terms  remain 

(^(m)  = 0j+(A2  . (A.15) 


) 

S-  l)a>  + 0)  + n(k-C-  1)] 

0)  ■ (A.13) 


where 


01  = E 


02  - E 


L-l 

2 [Wj{k)Wj(k  + m)l  A2sin(w(k-j-l)  + 0)sin(cj(k  + m-i-l)  + 0)J 
ij,=0 

(A.16) 
L-l 

^ Wj(k)Wj(k  + m)n(k-j- l)n(k  + m-i- 1)  | . (A.17) 

ij=0 


In  reference  36  it  is  shown  that  the  misadjustment  noise  V.(k)  is  highly  correlated 


in  time  and  consequently  (compare  (8))  that 

E(Vj(k)  Vj(k  + m))  = #if25(j_j) 


(A.18) 


A 


T 
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This  and  (7)  yield 


E[(Wj(k)  Wj(k  + m)]  = w*  w*  + 6(i  - j)  . 

We  now  calculate  0j.  Substituting  (A.  13)  and  (A.  19)  in  (A.  16),  we  have 


(A.19) 


^sin(c«;(k-j-  1)  + 0)  sin  (u)(k-j  + m-  1)  + d)  + E(r(k  1 6))  E(r(k  + m 1 0)) 
j=0 


Using  the  identities 

■ / j cos  b cos  (2a  + b) 
sin  a sin  (a  + b)  = — 

and  (A.  1 ) and  (A. 3)  we  obtain 

^ j = cos  c«jni(L  + terms  « L)  + E(r(k  | 6))  E(r(k  + m)  6)) 


(A.20) 


(A.21) 


^ rt^Lcoscom  + E(r(k|0))E(r(k  + m)0))  . 

Since  the  input  noise  n is  assumed  to  be  white,  (A.  17)  and  (A.19)  yield 
L-1 

02=  2 (wj^  W£  + 6(k  - C)/i$^)  6(m  - C + k) 
k,C=0 

L-l-|m| 

k=0  i 


|m|^  L 


From  (11)  and  (A.21)  we  have 


wj^  '^'k+lml  “ — 1 ) cos  cj(k  + m + 1 )) 
Ir 


(A.22) 


(A.23) 


(A.24) 


■ (cos  cum  + cos  u)(2k  + m + 2)) 


(A.25) 


Substituting  (A.25)  into  (A.24)  and  using  the  identity  (A.2)  we  obtain 
02  = 02  ) lm|  < L 


0 Mml  > L , 


(A.26) 
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where 


where 


Let 


Then  substitution  of  (A. 28)  in  (A.l  1 ) yields  four  terms,  each  containing  a double  sum  over 
6,k=0,  . . K-1  and  three  of  the  terms  restricted  to  |k-C|  < L-1.  Thus 


'>222? 


(A.34) 


where 


oj  - ^ a6(C  - k)  cos  cok  cos  oj2  , 
|C-k|<L 

K-1 

2 b cos  (K  - k)  cos  cuk  cos  coS  , 
£,k 


(A.35) 


(A.36) 


22 


^3  - ^ c (L  - I C - k I)  cos  co(C  - k)  cos  cok  cos  coE  , 
|E-k|<L 


~ 2/  ~ ■ 
|fi-k|<L 

If  we  make  the  substitution  n = £-k,  then 


M-l 

K-|n|-l 

2 " 

2 

2 . 

|E-k|<L 

n-l-M 

k=0 

K-l 

K-l 

2 - 

2 

2 • 

E,k 

nE=0 

k=0 

where  M = min(L,K)  as  in  (A.5). 

M-l 

K-|n|-i 

‘^l=a  2 

2 

6(n)  cos  cok  cos  co(n  + k) 

l-M 

k=0 

K-1 

= a (cos  cjk)^ 
k=0 


But  (cos  wk)2  = j + cos  2<ok  ^ 


2^  aK 


a,- 


We  now  calculate  o^: 


K-l  K-l 

^2  ~ co(C-k)  cos  coE  cos  (cok) 

E=0  k=0 


K-l  K-l 

“ 4 ^ ^ ^ "*’  2co(E  - k)] 


E=0  k=0 


Under  (A.3) 


K-l 

K-l 

K-l 

K-l 

2 

2 cos  2co(E  - k) 

<2 

7 J cos  co(E  - k) 

E=0 

E=0 

1 o 
>1 

1 cy 

k=0 

(A.37) 


(A.38) 


(A.39) 


« K • K = K^ 


(A.40) 
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Similar  expressions  for  the  remaining  terms  yield 


a?  = — + terms  « 

2 4 


(A.4I) 


To  calculate  a^' 

M-1  K-ln|-l 

*^3  “ c(L  - ln|)  cos  con  cos  cok  cos  co(k  + n) 

n=-M+l  k=0 

M-1  K-n-1  K-1 

= 2 ^ ^c{L  - n)coscon  coscok  cos  co(k  + n) -cL^cos2  cok  (A.42) 

n=0  k=0  k=0 

We  note  using  (A.  1 ) to  (A.3)  that 
K-1 

^ cos^  cok  « K 
k=0 

We  further  note  that 


M-1 

K-n-1 

K-n-1 

2 

z 

t 

Lcos2co(k+n)  < 2ML  /^cos  ixcojk  + n) 

n=0 

k=0 

k=0 

and  that 

M-1 

k-n-1 

2 

s 

2 

ncos  2co(k  + n) 

n=0 

k=0 

M K-n-1 

= 2 ^ ^ (n  cos  2cok  cos  2con  - n sin  2wk  sin  2con  «4m2k<4MKL  . 

n=0  k=0 

(A.44) 

Similar  expressions  hold  for  cos  2con  and  cos  2cok.  The  substitution  of  identity  (A.7b)  and 
(A. 42)  and  the  use  of  approximations  (A. 43),  (A.44),  etc.,  yield 


M-1  K-n-1 

3 = 2 2 2 T (L-n)+  J (terms  «MLK). 


'3  ^ ^ 4^  4 

n=0  k=0 


(A.45) 
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Also 


and 


M-1 

K-n-1 

M-l 

s 

n 

L 2 (K-n)  = 

n=0 

k=0 

n=0 

M-1 

K-n-1 

M-1 

2 

Z " 

= (K  - n)n  = 

n=0 

k=0 

n=0 

(A.46) 


(A.47) 


Combining  (A.45)  through  (A.47),  we  have 
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where 


K^L 

k3 

1 4 

12 

' L^K  _ 

l3 

. 4 

12 

, K>L  . 

\ 4 12 

Since  is  approximately  independent  of  o),  we  must  have 

2 


(A.53) 


CALCULATION  OF 

K-1  K-1 

Ouv  “ 2 2 ~ ■ (A.54) 

k=0  e=o 

Substitution  of  (A.32)  in  (A.50)  yields  four  terms 
®uv  = “i  +02+ “3  ■^“4 

similar  to  (A.35)  through  (A.38)  except  that  cos  cuk  is  replaced  by  sin  tuk.  The  first  term 
gives 

aj  = a:6()2  - k)  cos  coC  sin  cok 

|C-kl<L 


= — ^ sin  2cuk  , |aj|  « — K 


(A.55) 


The  next  two  terms  contain  the  expression 


cos  co(8  - k)  cos  cj£  sin  cuk  = — [sin  2co(8  - k)  + sin  2u;6  + sin  2a;kl 

4 

As  in  (A.40),  (A.43),  and  (A.44) 

I (sin  2w(£  - k)  « , 

C.k 

M-1  K-n-1 

(sin2co(k  + n)  «MLK  , 

n=0  k=0 


(A.56) 


M-1  K-n-1 

2 ^ n sin  2w(k  + n)  « 2MLK 

n=0  k=0 

Similar  expressions  are  valid  for  sin  2co2  and  sin  2cok.  Thus 

0(2  “ sin  a;(6  - k)  sin  cj6  sin  cok 

£,k 

= — (terms  « K^)  ; 

4 

M-1  K-|n|-l 

“3  “ c (L  - |n|)  cos  cun  cos  cu(n  + k)  sin  cuk 

n=-M+l  k=0 

= - (terms  « LMK)  ; 

4 

M-1  K-|n|-l 

0(4  = ^ ^ c d sin  cun  cos  cu(n  + k)  sin  cuk  . 

n=-M+l  k=0 


Once  again 


(A.57) 


(A.58) 


(A.59) 


(A.60) 


sin  cun  cos  cu(n  + k)  sin  cuk  = ^ [ 1 + cos  2cu(n  + k)  - cos  2cun  - cos  2cuk] 


\oi^\  < — (2MK  + terms  « MK) 


«-(LMK)  . 
4 


(A.61) 


Comparing  (A.55),  (A.58),  (A.59),  and  (A.61)  with  aj,  a^,  aj,  and  a\  we  see  that 
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